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The binding energy and the corresponding wave function of excitons in GaAs-based finite square 
quantum wells (QWs) are calculated by the direct numerical solution of the three-dimensional 
Schrodinger equation. The precise results for the lowest exciton state are obtained by the Hamilto¬ 
nian discretization using the high-order finite-difference scheme. The microscopic calculations are 
compared with the results obtained by the standard variational approach. The exciton binding 
energies found by two methods coincide within 0.1 meV for the wide range of QW widths. The ra¬ 
diative decay rate is calculated for QWs of various widths using the exciton wave functions obtained 
by direct and variational methods. The radiative decay rates are confronted with the experimental 
data measured for high-quality GaAs/AlGaAs and InGaAs/GaAs QW heterostructures grown by 
molecular beam epitaxy. The calculated and measured values are in good agreement, though slight 
differences with earlier calculations of the radiative decay rate are observed. 


I. INTRODUCTION 

Excitons in bulk semiconductors and heterostructures 
have been under intensive study for many years since 
their discovery in 1952 [T]. One of the important charac¬ 
teristics of an exciton is the binding energy caused by the 
Coulomb interaction of the electron and the hole |M1. 
Although in bulk semiconductors this energy is relatively 
small, typically lower than the lattice vibration energy 
at room temperature, in the semiconductor heterostruc¬ 
tures it can increase significantly, up to four times [2]. 
Together with the binding energy, the radiative proper¬ 
ties of an exciton are characterized by another important 
parameter, the radiative decay rate [4] or the oscillator 
strength f 7j, which is defined by the exciton-light cou¬ 
pling. Since the discovery of the giant exciton oscillator 
strength (the Rashba effect [8]), the exciton states and 
exciton-light coupling have been drawing much attention, 
see, for example, Refs. [2HEE]. Recent developments of 
the microcavities mm open up new frontiers for con¬ 
trolling the exciton-light coupling efficiency. 

Precise experimental determination of the exciton 
binding energy is a quite complicated problem [8]. The 
exciton transition can be shifted and broadened due to 
defects and roughness of the interfaces in a heterostruc¬ 
ture. The edge of optical transitions between uncou¬ 
pled electron and hole is usually inexplicit, therefore the 
spectroscopic determination of the exciton binding en¬ 
ergy leads to significant uncertainties. The most accu¬ 
rate methods seem to be a spectroscopic study of optical 
transitions to different exciton states (Is, 2s, ... ). The 
experimental study of exciton-light coupling and quanti¬ 
tative determination of its characteristics are also a non¬ 
trivial problem. The coupling gives rise to an energy shift 
and a radiative broadening of exciton lines |J] , which can 
be, in principle, measured using, e.g., the steady-state 
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reflectance spectroscopy, see Refs. BM5] and references 
therein. This approach, however, is applicable only for 
the high-quality heterostructures when the non-radiative 
processes do not dominate over the radiative one. An¬ 
other approach is the time-resolved spectroscopy using 
photoluminescence kinetics, pump-probe, or four-wave 
mixing methods for determination of the radiative decay 
time, see, e.g., Refs. miMIS!- Similar problems appear 
in this approach if the quality of the heterostructure is 
not high enough. 

The theoretical modeling of the exciton in a bulk semi¬ 
conductor is usually carried out in the framework of the 
hydrogen model [7j. In this model, the motion of an ex¬ 
citon as a whole and relative motion of the electron and 
the hole are separated that leads to the two independent 
Schrodinger equations (SEs) for the center-of-mass and 
intrinsic motions. The solution of the first equation is 
the plane wave, whereas the latter equation is similar 
to the SE for a hydrogen atom and can be analytically 
reduced to the one-dimensional problem. In this case, 
the exciton binding energy is easily calculated for var¬ 
ious semiconductors. For example, in the case of bulk 
GaAs semiconductor, this energy is R x = 4.26 meV [20} . 

The hydrogen model becomes unsuitable for semicon¬ 
ductors with degenerated valence band. In this case, 
the Kohn-Luttinger Hamiltonian [2T|, which couples the 
center-of-mass and intrinsic electron-hole motions, is 
more appropriate. As a result, the variables of the exci¬ 
ton SE cannot be separated and one has to consider the 
multi-dimensional SE [20]. A study of the exciton in a 
quantum well (QW) meets further complications of the 
problem. Even in the simplest case, when only the diag¬ 
onal part of the Kohn-Luttinger Hamiltonian is included 
in the problem, the presence of the QW potential requires 
one to consider at least the three-dimensional SE, which 
cannot be solved analytically. 

The exciton states in QW have been theoretically and 
numerically studied by many authors. The pioneer¬ 
ing work of Miller et al. [22|| has considered the two- 
dimensional exciton approximation for narrow QWs and 
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obtained the exciton binding energy close to 4R X . Fol¬ 
lowing works of Bastard et al. [2] and Greene et al. [23] 
employed the variational approach with various types of 
trial wave functions. In these works, the QWs with finite 
barriers as well as with infinitely high ones were studied. 
The binding energies have been calculated for different 
QW thicknesses and the breakup of the 4R X limit has 
been obtained for excitons in the QWs with finite bar¬ 
riers. Similar results have been obtained in the paper 
of Leavitt and Little [24] who have applied the method 
based on the adiabatic approach. It is worth noting 
that these authors have also studied the indirect exci¬ 
tons. More recent work of Gerlach et al. [25] describes 
a phenomenological model for the binding energy of an 
exciton in a QW. The authors compare the results of 
their model with the energies obtained by the variational 
approach with the prescribed trial function for the ex¬ 
citon wave function. The comparison shows a certain 
advantage in accuracy of the variational approach over 
the phenomenological model. 

The trial wave functions used in the variational calcu¬ 
lations of the exciton binding energy can be accurately 
chosen for two extreme cases. One of such cases is the 
wide QW, when the width, L , is much larger than the ex¬ 
citon Bohr radius, a#, so that the motion of the exciton 
as a whole can be separated from the intrinsic motion of 
the electron and the hole. In this case, the trial func¬ 
tion is a product of the hydrogen-like ls-function and 
a standing wave describing the quantization of the exci¬ 
ton motion in the wide QW. Another case is the narrow 
QW, in which the quantum confined energy of the elec¬ 
tron and the hole is much larger than their Coulomb cou¬ 
pling energy. Then, the separately obtained wave func¬ 
tions of the quantum confined electron and hole states 
are appropriate approximations [25]. The part of wave 
function taking into account the Coulomb interaction is 
usually modeled by a simple function with some parame¬ 
ters si m- So far, the precision of such approximations 
has not been studied in detail. 

The calculations of the radiative decay rate (or the os¬ 
cillator strength) have been carried out in several works, 
see, e.g., Refs. PS!F28] . The robust theoretical cal¬ 
culations have been done by Iotti and Andreani f25] 
and D’Andrea et al. |30j . where the general minimum 
of the oscillator strength for the GaAs/AlGaAs and In- 
GaAs/GaAs QWs has been found to be at QW width 
L ~ 2.5 as- This minimum defines the transition from 
the so-called weak confinement (L as) to the strong 
confinement (L < as), when the electron and the hole 
are separately localized. The experimental determina¬ 
tion of the radiative characteristics is widespread, see, 
e.g., Refs. ,10] I3H135] and references therein. Recent 
measurements of the radiative decay rate have been car¬ 
ried out by Poltavtsev et al ., see Refs. mug. In these 
papers, the general theoretical behavior has been exper¬ 
imentally confirmed, but the spread of the experimental 
measurements and the shortage of the available theoret¬ 
ical calculations motivated us to fulfill this deficiency. 


In the present paper, we provide the results of numer¬ 
ical solution of the SE for excitons in QWs with degen¬ 
erated valence band. Partial separation (over two vari¬ 
ables) of the center-of-mass motion and cylindrical sym¬ 
metry of the problem allowed us to reduce the initial SE 
to the three-dimensional one. We numerically solve the 
three-dimensional SE for the heavy-hole exciton in QWs 
of various widths and calculate the exciton binding en¬ 
ergy and the radiative decay rate. 

The numerical solution of the problem has been done 
for the GaAs/AlGaAs and InGaAs/GaAs QWs, which 
are widely experimentally and theoretically studied now 
as the model heterostructures [33, _3Bi. The obtained so¬ 
lutions are compared with results of the variational ap¬ 
proach with the trial function proposed in Ref. [25j . The 
results of the numerical experiments for different widths 
of the QWs are discussed in detail. 

We also experimentally measured the reflectance spec¬ 
tra for several high-quality heterostructures with In¬ 
GaAs/GaAs and GaAs/AlGaAs QWs grown by molec¬ 
ular beam epitaxy. An analysis of exciton resonances in 
the spectra using the theory described in Ref. [J allowed 
us to obtain the radiative decay rates for excitons in these 
structures and to compare them with the numerically ob¬ 
tained results. 

The paper is organized as follows. In Section 2, we 
derive the three-dimensional SE from the complete exci¬ 
ton Hamiltonian. The direct and variational methods of 
numerical solution are described in Section 3. The next 
section presents results on the exciton binding energy 
obtained by these two methods. Section 5 contains the 
numerically obtained results for the radiative decay rate, 
their comparison with experimental data, as well as an 
analysis of the decay rate for the narrow and wide QWs 
in the framework of simple models. The comparison of 
radiative decay rates and corresponding wave functions 
obtained by two numerical methods are discussed in Sec¬ 
tion 6. Section 7 summarizes main results of the paper. 
Additionally, some details of the numerical methods and 
their accuracies are given in the Appendix. 

II. MICROSCOPIC MODEL 

The exciton in a QW is described by the SE with 
Hamiltonian 

H = T e +T h - - t + I4(r e ) + Vft(r h ). (1) 

r h 

Here, indices e and h denote the electron and the hole, 
respectively. We introduce the kinetic operators for the 
electron, T e , and for the hole, T/, the relative electron- 
hole distance, |r e — r h |, the electron charge, e, and the 
dielectric constant, e. The square QW potential is given 
by 

V j (r J ) = V j e(^-^j, ( 2 ) 



3 


where j = e,h. Here 0(z) is the Heaviside step func¬ 
tion, L is the QW width, and Vj are the conduction- and 
valence-band offsets. 

In the effective-mass approximation m, the kinetic 
operator of the electron in the conduction band is explic¬ 
itly given as 


T e — Eg + 



(3) 


where E g is the energy band gap, k e = —ih\7 e is the 
electron momentum operator, and m e is the electron ef¬ 
fective mass. The kinetic term of the hole in the valence 
band is given by the Kohn-Luttinger Hamiltonian eh , 
which can be written as 


T h = — 
mo 


^7l + ~0~ + khyjy + kl z Jz) 


-273 (E k [ ln } { J rn , Jn\ 


(4) 


where 7 1 with l = 1,2,3 are the Luttinger parameters, 
khm, m = x,y,z, are the components of the hole mo¬ 
mentum operator, J m are the 4x4 angular momentum 
matrices [37:, {J m , J n } = (J m J„ + J„J m )/2 denotes the 
anticommutator, mo is the mass of the free electron. 

We consider only the diagonal part of the expres¬ 
sion Q, assuming that the nondiagonal terms contribute 
little to the energy states [20] ■ Therefore, we do not con¬ 
sider several effects extensively discussed in literature, 
see, e.g., Refs. (55, ESI- In particular, we ignore the 
heavy-hole - light-hole coupling. On the one hand, these 
effects give rise to a relatively small change of the car¬ 
rier and exciton energies and weak modification of the 
exciton-light coupling not exceeding experimental errors 
of the data obtained in our work. On the other hand, 
these simplifications allow us to get a precise numerical 
solution of the problem and to compare it to that ob¬ 
tained by the standard variational approach. So, up to 
the constant energy gap E g , the Hamiltonian |T]) has the 
form 


-Udiag 

+( 7 i T 272)/^ 


2 m e 


+ (71 ± 72) 


( k hx + k hy) 


2mo 


2mo e|r e - r h 
Introducing the effective masses 


+ V e (r e ) + V h (r h ). (5) 


yTlhxy 


m 0 

7 i ± 72 ’ 


m hz = 


m 0 

7 i T 272 ’ 


( 6 ) 


of the heavy hole (upper sign) and light hole (lower sign), 
respectively, we come to the Hamiltonian under consid¬ 


eration: 


| hz 

2 m hz 


-^diag — 
„2 


(khx + k hy ) 


e|r e - r h 


2 m e 2lJly lX y 

+ Re( r e) + Vh{ Fh)- 


(7) 


In our study, we pay attention only to the heavy-hole 
excitons. The theoretical analysis for the light-hole exci- 
tons is similar and differs only in Eq. (J 6 |. 

With the Hamiltonian |7|, we come to the standard 
six-dimensional SE, Hdiag'k = E^>, for the electron and 
the hole coupled by the Coulomb interaction JDj. The 
translational symmetry along the QW layer allows us to 
reduce this equation only to the four-dimensional one 
by separation of the center-of-mass motion in the (x, y)- 
plane. This motion is described by an analytical part 
of the complete wave function, 4/. The relative motion 
of the electron and the hole in the exciton is described 
by the part ip{x,y, z e , Zh) of the complete wavefunction, 
where x = Xh ~ x e , y = yh — y e - One more dimension is 
eliminated by taking advantage of the cylindrical symme¬ 
try of the problem and introducing the polar coordinates 
(p, cf> ) for description of the relative motion, since the 
Coulomb potential does not depend on </>. Representing 
the wave function in the form 


V>( x,y,z e ,z h ) 


j>{z e ,z hl p) 

P 


exp 


( 8 ) 


where k$ = 0,1,2,..., we proceed to the three- 

dimensional SE, which is numerically studied in the 
present paper. In Eq. ([8]), we introduce the factor 1/p 
in order to fulfill the cusp condition EH at p = 0. 

Since the light interacts mainly with the ground Is 
state of the exciton, we study the case when = 0. In 
this case, the SE under consideration is written as [2] 

( K - & == + V e (z e ) + V h (z h ) ]i>(z e ,z h ,p) 

\ ^VP +( Z e-Z h ) 2 J 

= E x ip(z e , z h , p) (9) 


where the kinetic term reads 


_ H 2 d 2 H 2 d 2 H 2 / d 2 Id 1\ 

2 m e dzl 2rrihz dz\ 2p \9p 2 pdp^ p 2 ) ’ 

( 10 ) 

and p = m e mhxy/ {fn e + nrihxy ) is the reduced mass in 
the (x, y)-plane. The obtained three-dimensional SE ([ 9 ]) 
cannot be solved analytically because the potential terms 
do not allow further separation of the variables. The 
approximate solution is also difficult for the finite QW 
potentials V e y, due to penetration of the wave function 
under the barriers. 

In our study, SE © is solved numerically and the en¬ 
ergy and corresponding wave function are obtained for 
QWs of various widths and compositions of the QW layer 
and barriers. Other parts of the complete wave function 
are analytically known, though they are not of the prac¬ 
tical importance. 
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TABLE I. Material parameters used for solving the eigenvalue problem ©• Energy gap mismatch denoted as A E g is calculated 
for AlGaAs heterostructure based on the data from Ref. 122 , for InGaAs, it is based on Ref. [30] . Right side of the table contains 
the effective masses for pure materials. Masses for the ternary alloys are obtained using the linear interpolation on x. Material 
parameters for AlGaAs and InGaAs heterostructures are listed here according to Ref. [25] and Ref. [44], respectively. 


Heterostructure 

X 

€ 

Ve/Vh 

A E g (meV) 


m e /mo 

mhz/mo 

ITlhxy / mo 

p/mo 

GaAs/ALrGai—a;As 

0.3 

12.53 

0.65/0.35 

365.5 

GaAs 

AlAs 

0.067 

0.150 

0.377 

0.478 

0.111 

0.242 

0.042 

0.093 

In^Gai-^As/GaAs 

0.02 

0.09 

12.53 

0.65/0.35 

0.55/0.45 

30 

139 

GaAs 

InAs 

0.067 

0.026 

0.350 

0.333 

0.111 

0.035 

0.042 

0.015 


III. NUMERICAL METHODS 


We performed the direct numerical solution of Eq. ^ 
for precise calculations of the exciton ground state en¬ 
ergy, E x . The exponential decrease of the wave func¬ 
tion at large values of variables allows us to impose the 
zero boundary conditions for the wave function at the 
boundary of some rectangular domain. The size of this 
domain varies from dozens QW widths (for small widths) 
down to several QW widths (for large widths). Therefore, 
the studied boundary value problem (BVP) is formed by 
Eq. § and the zero boundary conditions at p = 0, some 
large value of p and at large positive and negative values 
of the variables z e ^. Since SE |9]) is the standard three- 
dimensional partial differential equation of elliptic type, 
the direct numerical solution of the BVP is feasible us¬ 
ing available computational facilities. For this purpose, 
we employed the fourth-order finite-difference approxi¬ 
mation of the derivatives on the equidistant grids over 
three variables. The precise values of the studied quan¬ 
tities are obtained by the extrapolation of the results of 
calculations as the grid step goes to zero. The details on 
the numerical scheme and theoretical uncertainties are 
described in the Appendix. 

The nonzero solution of the homogeneous equation © 
with trivial boundary conditions can be obtained by 
the diagonalization of the matrix constructed from the 
Hamiltonian of Eq. ©• The obtained square matrix is 
large, nonsymmetric and sparse. The typical size of the 
matrix is of the order of 10 6 , so we keep in the calcula¬ 
tions only nontrivial matrix elements of a few diagonals. 
The diagonalization of such a matrix is difficult, however, 
a small part of the spectrum can be easily obtained. Us¬ 
ing the Arnoldi algorithm (45], we have calculated the 
lowest eigenvalue of the matrix and the corresponding 
eigenfunction. As a result, the ground state energy, E x , 
and the corresponding wave function have been obtained 
for various widths of QW. 

One of the alternative numerical methods is the ap¬ 
proximate solution of the three-dimensional SE © by 
the variational approach. This technique has been ap¬ 
plied by many authors, see Refs. mmmm- In the 
framework of the variational approach, the ground state 
energy of the system with the Hamiltonian H x is deter¬ 


mined by the minimization of the functional 

(ip\H x \ip) 


F = 


(V#) 


( 11 ) 


with respect to some free parameters of the trial wave 
function ip. The Hamiltonian H x corresponds to the left- 
hand side part of Eq. ([9]). For numerical calculation of the 
integrals in Eq. © one has to define the trial function. 
In Ref. I25j . the trial wave function having the form 

ip{z e ,z h ,p) = ip e (z e )iph(z h ) exp ( —— vV + A (z e - z h ) : 

\ as 

( 12 ) 

is applied. Here, the functions ip e (z e ) and iph(zh) are the 
ground state wave functions of the free electron and free 
hole, respectively, and a, A are the varying parameters. 

The shortcoming of the described variational approach 
is that the trial wave function has the prescribed form. 
Moreover, it assumes the partial separation of the vari¬ 
ables, whereas the Coulomb potential in Eq. ([9]) does not 
allow that separation. Of course, one can define even 
more complicated trial functions [231 [43] , but asymptoti¬ 
cally (for small and large widths of QW) they have to be 
reduced to the function (12). This fact as well as the nu¬ 


merical simplicity of this anzatz provoked many authors 
to use it for calculation of the exciton binding energy, 
R x . However, the accuracy of the results has not been 
studied in detail so far. The obtained wave function has 
also not been applied to calculate the exciton radiative 
characteristics. 

We apply the described numerical algorithms for solv¬ 
ing the SE with parameters for the GaAs/ARGai-^As 
and In^Gai-ajAs/GaAs QWs, which are widespread in 
the contemporary experimental studies. The parameters 
are general for such types of the QWs and, in particu¬ 
lar, simulate the typical ratio of the band offsets at the 
GaAs/Al 0 . 3 Gao. 7 As interface: V e /Vh ~ 2. In the calcula¬ 
tion, we have used the heavy-hole ARGai-^As parame¬ 
ters reported in Ref. [25] as well as the In^Gai-zAs ones 
in Ref. |44j . They are presented in Table [T] Parameter x 
denotes the concentration of A1 (In) in the Al(In)GaAs 
solid solutions, AE g (x) is the difference of the energy 
gap E g (x) for x = 0 and x > 0. For InGaAs, ratios of 
potential barriers were taken V e /\ 4, = 65/35 for concen¬ 
tration x = 0.02 and V e /Vh = 55/45 for x = 0.09. For 
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simplicity, we neglect the difference of the electron and 
hole masses in the QW and in the barrier layers. We 
ignore the small discontinuity of the dielectric constant 
at the QW interfaces as well H5] , 


IV. EXCITON BINDING ENERGY 


The exciton binding energy, R x . is defined by the exci- 
ton ground state energy, E x , with respect to the quantum 
confinement energy of the electron, E e , and the hole, Eh, 
in QW. This relation is given by the formula 

Rx = E e + Eh — E x , 

which characterizes the energy of the Coulomb interac¬ 
tion of the electron and the hole trapped in the QW po¬ 
tentials V e and 14, respectively. Energies E e and Eh 
are obtained from solution of the corresponding one¬ 
dimensional SEs for the electron and the hole in QW [3]. 

The direct microscopic calculations of the exciton bind¬ 
ing energy have been carried out for various widths of the 
GaAs/Alo 3 Gao. 7 As QW. The concentration of aluminum 
in GaAs has been chosen to be x = 0.3 in accordance to 
the concentration of the grown samples. Fig. [T] shows 
the exciton binding energy, R x , obtained from the mi¬ 
croscopic calculations and the variational approach as a 
function of the QW thickness. Surprisingly, the results 
of two methods coincide with the precision of 0.1 meV 
for QW widths L > 2.5 nm. This result means that the 
trial function (121 is a good approximation of the exci¬ 


ton wave function for the lowest state. The variational 
approach gives a bit smaller binding energy, which is the 
expectable result because the minimum of functional GD 
is achieved for exact exciton wave function rather than 
for the trial function. 

The uncertainties of the microscopic calculations of 
R x come from the calculated values of E x because the 
quantum confinement energies, E e and Eh, can be cal¬ 
culated with an arbitrary precision. We obtained the 
relative uncertainty of R x to be smaller than 1% for 
1.0 < L < 2.5 nm and decreasing down to values smaller 
than 0.1% for wider QWs. Although the uncertainty is 
much smaller than the difference of the results of about 
0.1 meV obtained by two numerical methods, this dis¬ 
crepancy seems to be negligible for comparison with con¬ 
temporary experimental data. Therefore, both methods 
can be successfully used for determination of the exciton 
binding energy. 

Interestingly enough, the maximum binding energy of 
9.71 meV is achieved by the direct microscopic solution at 
the QW width of about 3 nm. The variational approach 
gives for this QW E x = 9.66 meV at the optimal param¬ 
eters of the trial function (12) are a = 1.438, A = 0.3125. 


In contrast to the case of QWs with infinitely high barri¬ 
ers, for this case, the binding energy decreases for thinner 
QWs due to the penetration of the carriers into the barri¬ 
ers. In the other limit of very thick QWs, we obtained the 
values approaching the free exciton binding energy for the 



FIG. 1. Exciton binding energy, R x , obtained by the direct 
numerical solution of Eq. (|9| (blue curve) and by the varia¬ 
tional approach (red curve) for GaAs/Alo. 3 Gao. 7 As QWs as 
a function of the QW width. 


bulk crystal. From Fig. [TJ one can see that, for our case, 
this energy is R x = 4.10 ± 0.03 meV that is slightly less 
than reported in Ref. [2D.;. For the variational approach, 
the parameters of the trial function (12) for the QW of 
width L = 200 nm, where the exciton can be treated as 
the free one, are a = 1.0604, A = 1.1858. These val¬ 
ues are close to the expectable ones for the bulk crystal: 
CToo — I; Aoo — I SI- The overall behavior of the binding 
energy for QW widths 1 < L < 150 nm can be approxi¬ 
mated with the accuracy of 0.07 meV by simple formula 
R X {L) = (3.5431L 2 +146.694L)/(L 2 + 11.1758L+6.0473). 


Together with the exciton ground state energy, in the 
direct microscopic calculation we have obtained the cor¬ 
responding wave function. In the wide QWs, the wave 
function ip(z e , Zh, p) can be presented as a function of 
relative coordinate, z = z e — Zh, and center-of-mass co¬ 
ordinate, Z = ( m e z e + nihZh) /(m e + m^). The slices of 
\ip(Z, z, p)\ 2 as functions of £ and p for three different 
center-of-mass coordinates Z are shown in Fig. |2j[a) for 
the QW width L = 150 nm. These slices show the proba¬ 
bility distribution for relative distance between the elec¬ 
tron and the hole in the exciton. The center plot shows 
the exciton localized in the center of the QW while the 
side plots present the exciton near the QW interfaces. 
The magnitude of | tl>(Z, z, p)| 2 for the side plots is in sev¬ 
eral orders smaller than for the center plot. 

As seen from Fig. [2j[a) , the probability distribution at 
Z = 0 nm is spherically symmetric and reveals some dis¬ 
tortion near the QW interfaces. The exciton in a wide 
QW can be considered as slowly propagating across the 
QW that is its kinetic energy is small compared to the 
Coulomb energy. Therefore an adiabatic approximation 
can be used, when the electron-hole relative motion in the 
exciton is considered to be much faster than the center- 
of-mass motion. In the framework of this approximation, 
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the obtained distortion can be treated as an appearance 
of a static dipole moment of the exciton when it encoun¬ 
ters the barrier. 


z (nm) 


z (nm) 


z (nm) 



Z (nm) 


decreases. 


V. EXCITON-LIGHT COUPLING 


A. Radiative decay rate 


Exciton-light coupling is usually characterized by ei¬ 
ther the radiative decay rate or the oscillator strength [5J 
0 m izzi- The radiative decay rate, To, characterizes 
the decay of electromagnetic field emitted by the exci¬ 
ton after the pulsed excitation: E[t) = if(0) exp(—Tot). 
A consistent exciton-light coupling theory is presented, 
e.g., in the monograph of Ivchenko [J|. It provides the 
folowing expression for Tq: 


27 iq / e\p cv \ 
He \m 0 Wo 


/ *( 2 )e*p (i , 2 )<fc , 

— oo 


(13) 



FIG. 2. (a) The slices of the exciton wave function squared, 

\ip{Z, z, p)| 2 , as functions of relative coordinates ( z,p ) for 
three different center-of-mass coordinates Z = —69, 0, 69 nm 
of the exciton. The QW width, L, equals to 150 nm. The co¬ 
ordinates of the center-of-masses of the excitons are depicted 
by small crosses. The QW interfaces are drawn by the dashed 
vertical lines. The probability for the electron and the hole 
to be at distance r = y/ z 2 + p 2 is shown by the false colors 
in the logarithmic scale, (b) The exciton dipole moment as a 
function of the center-of-mass coordinate. 

The exciton static dipole moment has been extensively 
studied for indirect excitons in coupled QWs, see, e.g., 
Refs. g7HH3 and references therein. In our case of the 
direct exciton, it has been calculated using the standard 
definition: D x = e (il)\z\4>). The dependence of the dipole 
moment on the exciton position in the QW is presented 
in Fig. [2] (b). It is clearly seen that the maximum dipole 
moment is achieved near the barrier of the QW, where the 
wave function undergoes the maximum distortion. When 
the wave function is spherically symmetrical (Z = 0), the 
dipole moment is zero and increases in absolute value as 
the exciton approaches the barriers. If the exciton is very 
close to the barrier, the wave function penetrates under 
the barrier, distortion diminishes and the dipole moment 


where q = y/euj/c is the light wave vector, ojq is the exci¬ 
ton frequency, \p cv \ = (u v \e ■ p|u c ) is the matrix element 
of the momentum operator between the single-electron 
conduction- and valence-band states, and ^(z) = ip(z = 
z e = z h ,p = 0). 


The simplification of Eq. (13 1 used in Refs. [23 125j 
takes into account that the QW width is much smaller 
than the light-wave length 2ir/q. It allows one to replace 
exp(iqz) in the integral by unity. In this case, the ra¬ 
diative decay rate T 0 is closely related to the oscillator 
strength per unit area [27], f/S, by the formula 


r n = 


/ 

nrriQC \S 


(14) 


where n is the refraction coefficient (y/e = n + ik). 

The wave function obtained from the microscopic cal¬ 
culation allowed us to calculate the radiative decay rate 
T 0 of the exciton ground state according to Eq. (131. 


We calculated T 0 for GaAs/Al 0 . 3 Gao. 7 As QWs of var¬ 
ious widths from 1 nm to 300 nm. In the calcula¬ 
tions, \p cv \ 2 = m 0 -E p /2, where E p = 28.8 eV for GaAs 
and E p = 21.5 eV for In As are taken from Ref. |44j . 
The exciton frequency wo is calculated using bandgap 
Eg = 1.520 eV for GaAs [4] and parameters listed in 
Tab. [U 

Fig- H shows the radiative decay rate in energy units, 
ATo, as a function of the QW width. The radiative decay 
rate reaches its maximum at the QW width of about 
130 nm, that approximately corresponds to the half of 
the light wavelength in the QW material, A(GaAs) = 
A(vac)/n(GaAs) = 230 nm, where n(GaAs) = 3.6 is the 
refractive coefficient of GaAs at the photon energy huj = 
E g ( GaAs). So, this maximum of T 0 corresponds to the 
maximal overlap of the exciton wave function <£>( 2 :) and 
the light wave, see Eq. (131. 


As the QW width decreases, T 0 also decreases due to 


the diminish of the overlap integral in Eq. (13). For 
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small QW widths, however, To grows that correlates with 
the increase of the exciton binding energy (compare with 
Fig. 0 . Therefore, we may suppose that this maximum 
of To is caused by squeezing of the exciton in the narrow 
QWs by the QW potential. The squeezing gives rise to 
increase of the probability to find electron and the hole 
in the same position (z e = Zf, and p = 0). Between these 
maxima of To there is a minimum for the QW width 
of about 30 nm, which corresponds to the exciton Bohr 
diameter. The presence of such minimum was earlier 
pointed out by Iotti and Andreani [29]. 

We compared our results with two simple approx¬ 
imations: the exciton in the bulk semiconductor for 
wide QWs and two-dimensional approximation for nar¬ 
row QWs. Both the approximations are shown in Fig. [3] 
by various curves. The wave function of the exciton in 
the bulk semiconductor is given as [4] 


ip{Z,r) = 


\Jna%L* 


COS (l * Z ) 6XP 


a-B ) 


(15) 


where r is the electron-hole distance, L* = L — 2 Lj is the 
effective QW width obtained from the QW width L and 
the dead layer Ld [45] ;3U| |Mj. The part of the wave func¬ 
tion ( p~5| ) depending on Z is represented only by function 
cos ( 7 tZ/L*). It leads us to the conventional definition 
of the dead layer, which is the distance from the QW 
interface to the point where the cosine approximation 
of the exciton wave function becomes zero. This defini¬ 
tion of the dead layer is illustrated by the inset of Fig. [3] 
(b). Substituting the wave function (15) in Eq. (13), one 
can calculate the radiative decay rate in a QW with in¬ 
finite barriers and constant dead layer Ld = 14.6 nm. 
It is shown by dashed curve in Fig. [ 3 ] (a). One can see 
that this approximation is acceptable for the QW widths 
L > 140 nm. 

A better approximation of the dependence r 0 (A) can 
be achieved using the variable dead layer [43]. We ex¬ 


tracted the variable dead layer by fitting the function (15) 


at r = 0 to the numerically obtained 4>(,c) [see inset 
in Fig. [3] (b)]. Extracted values of Ld are shown in 
Fig- 1(b) by solid points. Dependence of Ld on the QW 
width L can be well approximated by the phenomeno¬ 
logical formula: Ld = a(l — exp (— L/L 0 )) + b, where 
a = 20.6±0.5 nm, L 0 = 70±2 nm, and b = — 5.7±0.7 nm. 
Using this approach and as = 14.55 nm, we obtained 
more accurate approximation of HTq{L) which is shown 
in Fig. 1 (a) by red solid curve. This approximation is 
appropriate for the QW widths down to 100 nm. 

For narrower QWs, both the bulk exciton wave func¬ 
tion and the idea of variable dead layer are no longer 
applicable. Instead, for thin QWs, we implemented the 
two-dimensional (2D) exciton approximation. The wave 
function of the 2D exciton has the form [4] 


ip(z e ,z h ,p) 


i>e{z e )^ z (z h ) 




(16) 


where p e s is the effective 2D exciton radius, ifj e (z e ) and 



QW width (nm) 



FIG. 3. (a) Radiative decay rate in energy units, /To, versus 
QW width, L. Blue dots are obtained in the precise micro¬ 
scopic calculation. Dashed line corresponds to To for the bulk 
exciton in a QW with infinite barriers and constant dead layer 
Ld = 14.6 nm. Solid red line is the calculation with the vari¬ 
able dead layer shown in (b). The green solid line corresponds 
to To in the 2D exciton model with the varying effective radius 
p e ff shown in (c). (b) The variable dead layer as a function of 
QW width. The inset: function $( 2 ) and its approximation 
by cos(-7 iZ/L*) for QW with L = 150 nm. The two-side ar¬ 
row marked 11 Ld” illustrates definition of the dead layer, (c) 
Effective 2D exciton radius, p e ff, extracted form the calcu¬ 
lations (solid points) for the QW widths 0-30 nm and its 
approximation by function ( |18| ) (solid curve) with parame¬ 
ters: Ci = 9.5 ± 0.4 nm, Li = 14 ± 1 nm, C 2 = 6.0 ± 0.5 nm, 


Z /2 = 2.4±0.3 nm. Das hed curves show contributions f aqu (L) 
and /pen(T) in Eq. Jl8l. 


iph{zh) are wave functions of the free electron and free 
hole in a QW with finite barriers. Function (161 is the 
solution of eigenvalue problem (|9|) with 2D Coulomb po¬ 
tential, — e 2 /(ep). In the true 2D exciton problem, p e g is 
the 2D Bohr radius 


P2D 


h 2 e 
2 /re 2 ’ 


(17) 


which is twice smaller than as- 

The exciton in the QWs with finite heights of the bar- 





















riers does not reach two-dimensional limit due to pen¬ 
etration of the exciton wave function into the barriers. 
Therefore, we consider p e g as a characteristic parame¬ 
ter. We fitted the wave function (16) to the numerically 
obtained wave function <I>(z) by varying p c g. Extracted 
effective 2D exciton radius values are presented in Fig. [3] 
(c). It is seen that p 0 Q decreases as L —► 4 nm and then 
increases again as L —> 0. This behavior is well described 
by a phenomenological function 


f(L) — fsqu(L) + fpen(L) 


P2D+C 1 (l-e- i / il ) 


+ 


C 2 e 


-l/l 2 


.(18) 


The first part of this function, / squ (L), decreases with 
L — > 0 down to the 2D-limit given by Eq. © > P'2B — 
7.9 nm for GaAs, and reflects the squeezing of exciton in 
narrow QWs. Another one, / pen (T), increases as L —» 0 
and reflects the penetration of exciton into the barrier. 
Using this dependency, we calculated HT o for the 2D exci¬ 
ton model, see Fig. [3] (a). As it is seen, this model, even 
with varying p e g, is adequate in terms of the radiative 
decay rate only for narrow QWs L < 15 nm. 

The described approximations are applicable for nar¬ 
row or wide QWs separately. The excitons in QWs of 
intermediate widths are not described by these models. 
Only the direct microscopic calculation provides the pre¬ 
cise values of the radiative decay rate for the wide range 
of QW widths. 


B. Experimental determination of To 


Exciton reflectance spectra can be used to obtain ra¬ 
diative decay rate. According to the theory summarized 
in Ref. [3], the amplitude reflectance coefficient of the 
QW with exciton resonance is given as 


«r 0 

r QW = - - 7r . p V 

UJq — oj — z(T + T 0 ) 


(19) 


where u >o is the renormalized exciton resonance frequency 
and u> is the frequency of the incident light. The nonra- 
diative decay rate T in Eq. (191 takes into account an 
additional broadening of the resonances due to nonra- 
diative processes. Reflectance tqw is strictly related to 
the reflectance coefficient of the whole heterostructure, 
R, which, in turn, can be measured in experiment. 

For the single QW heterostructure the relation between 
R and tqw can be found in Ref. 3]. In a structure with 
several QWs as those used in our study, the reflectance 
coefficient is generalized to 


r + Ej rqwj-e^ 

1 + r E , rQW jC' 0 ' 


( 20 ) 


Here r is the Fresnel reflectance coefficient from the sur¬ 
face of a heterostructure and <j> is the phase shift of the 
light wave reflected by the QW with respect to that re¬ 
flected by the structure surface. Effectively, Eq. (20) is 



Energy (meV) 



FIG. 4. (a) The experimental reflectance spectrum (blue 

curve) of the heterostructure with GaAs/AlGaAs QWs as well 
as the fit (red curve) of the spectrum. The ground heavy- 
and light-hole exciton states are denoted by X^h and Xih, 
respectively, (b) The same is for the heterostructure with 
InGaAs/GaAs QWs. 


the direct relation between the reflectance R and the ra¬ 
diative decay rate To- 

For comparison of the theoretical modeling with ex¬ 
perimental results, several high-quality heterostructures 
with QWs grown by molecular beam epitaxy have been 
selected for study of their reflectance spectra. The spec¬ 
tra were measured using a simple setup consisting of 
a white light source, a cryostat, and a spectrometer 
equipped with a CCD camera. Special precautions have 
been taken to accurately calibrate the absolute value of 
reflectance. For this purpose, a monochromatic light of 
a continuous-wave titanium-sapphire laser was used to 
measure the reflectance at a spectral point beyond the 
exciton resonances. 


Fig. ^demonstrates examples of the reflectance spectra 
(blue curves) for the GaAs/AlGaAs and InGaAs/GaAs 
heterostructures. We used the Eq. (20) to fit the spec¬ 
tra with To, r, wo and <j> as the fitting parameters for 
each resonance (red curves). The radiative decay rates 
in energy units, HT 0 , obtained for series of QW het¬ 
erostructures are presented in Tab. |TT] We fitted only the 
ground heavy-hole and light-hole exciton states. In the 
InGaAs spectra, we were able to reliably determine only 
the heavy-hole exciton resonances because the light-hole 
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ones are significantly detached by the strain |52| . The 


TABLE II. Radiative decay rate ftTo in energy units (/reV) 
together with the fit standard deviations. The data were ex¬ 
tracted from the measured spectra of samples with QWs of 
various widths. 


Structure 

QW width (nm) 

ATo of Xhh 

fiT 0 of X ih 

GaAs/AlGaAs 

14 

36.6 ±0.6 

14.9 ±0.7 


20 

35.4 ±0.4 

10.9 ±0.4 

InGaAs/GaAs 

2 

30.0 ±0.2 



3 

27.3 ±0.5 



3 

25.8 ±0.6 



30 

38.2 ±0.8 



35 

35.1 ±0.6 



40 

36.8 ±0.4 



95 

58.0 ±0.5 



shape of exciton peculiarities is known to be defined by 
phase 0j in Eq. (20 1 determined by the QW-to-surface 
distance in real sample. Relatively small peculiarities of 
the spectrum observed in Fig. |4|b) are, most probably, 
the excited quantum confined exciton states, which are 
beyond the scope of the present paper. 

Experimentally obtained radiative decay rates for 
heavy-hole excitons listed in Tab. [Tl] are compared to 
the calculated data in Fig. [5] In the calculations, we 
simulated GaAs/Alo. 3 Gao. 7 As and Ino. 02 Gao. 9 sAs/GaAs 
heterostructures with 365.5 meV and 30 meV band off¬ 
sets (see Tab. [I]), respectively. The calculations indi¬ 
cate that the radiative decay rate for GaAs/AlGaAs het¬ 
erostructure grows as L diminishes in the range L = 
5 — 20 nm. The experimentally studied GaAs/AlGaAs 
structures exhibit perfect agreement with the calculated 
results. In particular, they support the tendency of in¬ 
crease of AT 0 with decrease of the QW width in the range 
L = 5 — 20 nm. This tendency is also confirmed by many 
experimental data in Ref. El- 

Experimental results for InGaAs/GaAs heterostruc¬ 
tures demonstrate a trend of ATo to diminish with de¬ 
crease the QW width down to L sa 5 nm, which is in 
agreement with the calculations. At the same time, the 
experimental data are more spread around values ex¬ 
pected from computation. Large spread of the data is 
also observed in Ref. [15] . We explain this spread mainly 
by the indium concentration variation from one sample 
to another in the range 2 - 5%. We should note that 
high mobility of indium atoms during the growth pro¬ 
cess makes InGaAs heterostructures less predictable as 
compared to AlGaAs ones, in particular, due to the seg¬ 
regation effect [52] , 


C. AT 0 in shallow QWs 

The large difference in behavior of the radiative decay 
rate for GaAs/AlGaAs and InGaAs/GaAs in the range 



FIG. 5. Radiative decay rate for In^Gai-xAs/GaAs 
QWs with 2% and 9% of In (blue and red curves, respec¬ 
tively) as well as for GaAs/AlGaAs QWs (green curve). The 
empty squares and circles show the experimental data for In¬ 
GaAs/GaAs and GaAs/AlGaAs QWs, respectively. 


of small QW widths requires a particular analysis. We 
believe that this difference is related to different depth 
of the QWs. To check this assumption, we carried out 
computations of ATo for different concentrations x in the 
Dq.Gai-^As/GaAs heterostructures that results in dif¬ 
ferent heights of the QW barriers. In particular, for 
x = 0.02 and 0.09 the barriers are V e = 19.5 meV and 
V e = 76.5 meV, respectively. 

In the inset of Fig. [5j the radiative decay rates for 
IiUcGai-a-As/GaAs heterostructures with x = 0.02 and 
x = 0.09 and for GaAs/Al 0 . 3 Ga 0 . 7 As heterostructures 
are shown. It demonstrates the evolution of the radiative 
decay rate in narrow QWs with the growth of the barrier 
height. As the height increases with growing of x, the 
peak of HT 0 becomes more pronounced and its maximum 
is shifted to the lower QW widths. As we noted above, 
the peak at L « 5 nm for the GaAs/AlGaAs QWs is 
formed due to exciton squeezing by the QW potential. 

In the case of In^Gai-^As/GaAs QWs with small con¬ 
centration of indium, x = 0 . 02 , there is no peak of HTq at 
the small QW widths. This is an indication of the weak 
exciton squeezing due to the small QW potential depth. 
For the In^.Gai-^As/GaAs QWs with x = 0.09, some in¬ 
termediate behavior is observed due to the intermediate 
hight of the barriers. To understand this behavior, we 
applied the 2D exciton model described in Sect. |V A] and 
extracted p e g for In x Ga.i _ x As/GaAs QWs with 2% and 
9% of indium [see Fig. [ 6 ](b)]. The dependence p e g{L) for 
the shallow QW {x = 0.02) reveals the weak minimum, 
which is an indication of the weak squeezing of the ex¬ 
citon. Correspondingly, no maximum of ATo should be 
observed for these QWs. Indeed, the calculated width 
dependence of Ar 0 shows monotonic rise with L —> 0, see 
Fig.^a). Curve p e ff(L) for deeper QW {x = 0.09) shows 
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QW width (nm) 



QW width (nm) 


FIG. 6. (a) The calculated radiative decay rate the 

IiizGai-a;As/GaAs QWs with for 2% and 9% of indium 
(curves with points) in comparison with the 2D exciton model 
(smooth curves). Blue curves represent 2% of indium in the 
heterostructures, while 9% of indium are shown in carmine, 
(b) The fitted effective 2D exciton radius p e g (solid points) 
for 2% and 9% of indium as well as the phenomenological 
approximations (see the text). 


the noticeable minimum at L « 5 nm that points out to 
an exciton squeezing. As a result, a maximum of HTq 
appears [Fig. |6](a)]. 

The rapid growth of ATo at L —► 0 in the InGaAs/GaAs 
QWs with 2% of indium [Fig. |6][a)] is explained by the 
stronger penetration of exciton into the barriers as com¬ 
pared to the deeper InGaAs QWs with 9% of indium. 
Due to the penetration, the overlap of the exciton wave 
function <I>(z) and the light wave [see Eq. ©] rises and 
the radiative decay rate increases as L —► 0. 

In conclusion of this section, we compare our results 
of ATo for InGaAs/GaAs QWs with those of Ref. [50] for 
the same QW widths and indium concentration x = 0.09. 
The latter results have been also recently compared with 
experiment in Ref. m- Our calculations show very simi¬ 
lar dependence of the radiative decay rate in energy units, 
which, however, is shifted lower by 4-6 /ieV. Since the 
controlled precision of our calculations for L < 50 nm is 
better than 1.5 /ueV, our results are improved in precision 
with respect to the earlier calculations. 


VI. DISCUSSION 


In the previous sections, we have shown that the exci¬ 
ton ground state as well as the radiative decay rate To can 
be calculated by the direct numerical solution of the ex¬ 
citon SE. The binding energy is obtained quite precisely 
by the direct numerical solution as well as by the varia¬ 
tional approach with prescribed trial function (12). The 
good agreement of the results of two numerical methods 
indicates that the chosen trial function is appropriate for 
the wide range of QW thicknesses. It seems that, for each 
QW width, the varying parameters allow one to properly 
scale the trial function and, thus, to accurately simu¬ 
late the exact wave function of the exciton ground state. 
Therefore, besides the exciton binding energies, the wave 
functions obtained by two methods should also be very 
similar. 



FIG. 7. The radiative decay rate in energy units obtained for 
GaAs/AlGaAs QWs by the precise microscopic calculation in 
comparison to the results from the variational approach as a 
function of the QW width. 


In this context, it is interesting to compare the radia¬ 
tive decay rates obtained using these two methods. The 
exciton ground state wave function obtained by the varia¬ 
tional approach for the GaAs/Al 0 . 3 Gao. 7 As heterostruc¬ 
ture has been used to calculate r 0 and comparison of 
Ar 0 is presented in Fig. [7] One can see that the val¬ 
ues of ATo calculated by both methods differ slightly for 
the wide range of QW widths. This difference of about 
10%, however, is larger than the difference of the exciton 
binding energies at these QW widths. In particular, for 
L ~ 5 nm the variational approach gives greater radia¬ 
tive decay rate than the direct numerical solution. Nev¬ 
ertheless, such difference is generally not so important for 
comparison with experimental data because the spread of 
our experimental data as well as of recent measurements 
(see Refs. (TUCCSI) substantially exceeds this discrepancy. 
As a result, the variational approach with the proposed 
in Ref. [25] anzatz (12) is appropriate for calculation of 


















11 


the exciton binding energy as well as the radiative decay 
rate. 

Until now, we have considered the exciton binding en¬ 
ergy and radiative decay rate, which both are the inte¬ 
gral characteristics of the exciton states. Although the 
comparison showed us that these characteristics are sim¬ 
ilar, the wave functions obtained by different methods 
may differ in some details. This difference may affect 
other properties of excitons, e.g., sensitivity to a mag¬ 
netic field H]. So, for the future prospects we have com¬ 
pared the wave functions themselves. Fig. [8] shows the 
difference (in percents) of the exciton wave functions ob¬ 
tained by two methods as a function of p and z for several 
widths of the GaAs/AlGaAs QWs. It is clearly seen that 
the main difference of the wave functions is observed near 
the z = 0, p = 0 point. This difference for small L does 
not exceed 5% and rapidly decreases with the z and p 
rise. For large L , the difference is larger but also rapidly 
decreases as z and p grow. 



FIG. 8. The difference (in percents) of the wave functions 
of excitons in GaAs/AlGaAs QWs obtained by the precise 
microscopic calculation and by the variational approach as a 
function of two variables: p, z = z e = Zh for QW widths 
L = 5,20,40,100 nm (from bottom to top). The functions 
are shifted vertically for visibility. 


VII. CONCLUSIONS 

In the present paper, we have numerically solved the 
three-dimensional SE and obtained the exciton binding 
energy as well as the ground state wave function. The 
studied SE was deduced from the electron-hole Hamilto¬ 
nian with the Kohn-Luttinger term for the valence band 
taking into account only the heavy-hole excitons. The 
SE was solved by two methods: the direct one and the 
variational one. The variational method has been used 
by number of author who were studied this problem. The 


direct microscopic solution using the fourth-order finite- 
difference scheme has been carried out for the first time. 
It allowed us to precisely calculate the ground state and 
the exciton radiative decay rate. 

We obtained the radiative decay rates for 
GaAs/AlGaAs and InGaAs/GaAs heterostructures 
for various QW widths 1 < L < 200 nm and compound 
concentrations. We found that, since the different 
concentrations lead to the different magnitudes of band 
offsets, the behavior of the radiative decay rate for 
narrow QW widths (2 < L < 10 nm) strongly depends 
on the concentration. Increase of the concentration 
leads to the higher barriers and to growth of the peak of 
the radiative decay. For very narrow QWs (L < 2 nm) 
the radiative decay rate also behaves differently. If the 
QW is shallow, then the radiative decay rate grows as 
L -» 0. Instead, if the QW is deep, then it seems to be 
decreasing. 

The obtained wave functions were used to test some 
simple models, namely the 3D exciton model (exciton 
in the bulk semiconductor) and the 2D exciton model. 
We have shown that these models are applicable only for 
wide and narrow QWs, respectively. The parameters of 
the models and the limits of applicability were estimated. 

The comparison of the exciton binding energies ob¬ 
tained by the direct microscopic calculation and by the 
variational approach showed a very good agreement with 
the discrepancy of less than 0.1 meV for a wide range 
of QW widths. The analogous comparison of the radia¬ 
tive decay rates gave the similar overall behavior of these 
quantities, but larger differences (up to 10%). Neverthe¬ 
less, even for narrow QWs, this difference can be consid¬ 
ered as insignificant because the available experimental 
data are more spread for these QWs. As a result, both 
numerical methods can be successfully used for calcula¬ 
tion of the exciton binding energy and the radiative decay 
rate. 

The experimental measurements of the radiative de¬ 
cay rate (HWHM) for several QW widths presented in 
the paper are consistent with the results of the precise 
calculations. The comparison with the earlier calcula¬ 
tions showed slight systematic shift (of about 4 — 6 peV) 
of the earlier results with respect to our numerical data. 
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Appendix: Numerical details 


In our study, we employ the well known finite- 
difference (FD) approximation [55H55] of the deriva¬ 
tives in Eq. Q due to its robustness as well as sparse 
and relatively simple form of the obtained matrix equa¬ 
tion. We consider the BVP for Eq. © at the domain 


[—Z e / 2,Z e /2] x [— Zh/2, Z/j/2] x [0,1?] over variables z e , 
Zh : and p, respecitvely, and specify the homogeneous 
boundary conditions at the boundaries. The equidistant 
grids over each variable are introduced by the formulas 
= ^e/2, Zh,l. = ^A 2?i ^h/^i Pm “ ^Ap, 

where A Zeh = Z e>h /(N Ze ^ + 1), A p = R/{N P + 1) are 
the grid steps and the indices k. I , m go from 1 to some 
integer values IV 2e , N Zh , and N p , respectively. 


We use the central fourth-order FD formula for approximation of the second partial derivative of ip(z e , Zh, p) with 
respect to x e ; 

-tpk-2,1 ,m + 16 ip k ~ 1 ,l,m 30 'lj>k,l,m 4“ (-1 ,Z,m ^ 1pk+2,1 ,m 


12A? 


(A.l) 


Here, the unknown wave function on the grid ip(z et k, Pm) is denoted as ipk,i,m- The same FD formula is employed 
for the second derivative with respect to Zh- The wave function ip at the knots beyond the considered domain over z e 
and Zh is taken to be negligible due to its exponential decrease in this region. We apply the noncentral fourth-order 
FD formulas for approximation of the first and second partial derivatives of ip(z e , Zh, p) with respect to p\ 


3lpk,l,m 10t^fc,Z,m-}-l T 6lpk,l,m+3 T 1pk,l,m+ 4 

12A p ’ 

10 Ipk ,Z,m 15-zfa,; ,m+1 4 + Ipk,l ,m+5 

12Af 


(A.2) 


It allows us to satisfy the trivial boundary condition at 
p = 0 and to avoid knots p m < 0. At the knots p m > R , 
the wavefunction is also assumed to be zero. It should 
be noted that this assumtion is possible only if the con¬ 
sidered domain is large enough. 

In the calculations, the grid steps over each variable 
have been taken to be the same, A = A 2e = A Zh = 
A p , and multiply associated with the QW width. The 


formulas ( A.1|A.2 ) define the theoretical uncertainty of 
the numerical solution of order of A 4 as A A 0. However, 
the discontinuity of the square potential <[2j) at the QW 
interfaces decreases the convergence rate of the solution 
over z e and Zh to order of A 2 , whereas the convergence 
rate over variable p is kept ~ A 4 . 

In order to choose the appropriate numerical scheme, 
we compared the exciton ground state energy, E x , cal- 
culat ed using the second-order FDs and fourth-order 
FDs ( A.1|A.2 1 in Eq. ([9]). The convergence of the ex¬ 
citon ground state energy as a function of the square 
of the grid step, A 2 , is presented in Fig. [ 9 ] for AlGaAs 
and QW widths L = 5,100 nm. Although both schemes 
provide convergence to almost the same values of energy 
as A -A 0, the rate of convergence is different. For the 
energy obtained using the second-order FD, the rate of 
convergence considerably depends on the width of QW. 
For QWs of small widths (of order of the Bohr radius as), 
the linear convergence (with respect to A 2 ) is observed 
for both schemes, whereas for wide QWs the second-order 


scheme gives nonlinear convergence rate as A A 0. The 
energy obtained using the fourth-order FD shows the lin¬ 
ear dependence on A 2 for the whole range of the studied 
QW widths. This linear behavior is not changed with 
increase of the QW width. As a result, we performed a 
least square fit of the calculated energy and extrapolated 
it to A = 0. The uncertainties estimated from the least 
square fit are quite small. Thus, we obtained the pre¬ 
cise exciton ground state energy. The typical grid step 
reached in the calculations with the QW widths compa¬ 
rable with as is A = 0.25 nm whereas for wide QW the 
achieved grid step A = 0.5 nm. 

For the radiative decay rate in energy units, fiFo, calcu¬ 
lated using the ground state wave function obtained from 
the direct numerical solution, the convergence is not so 
perfect. The examples of the convergence rate for direct 
calculation of ATo are shown in Fig. [To] For different QW 
widths, the convergence rate is different. We were able to 
fit the convergence rate as A A 0 by the power function 
of A and estimate the uncertainty as the discrepancy of 
the calculated and extrapolated values. This discrepancy 
was smaller for narrow QW than for the wide QW. As a 
result, the uncertainty of our extrapolation for obtaining 
the precise radiative decay rate was estimated as 1.5 peV 
for QW width L < 50 nm and 3 peV for wider QWs. For 
the variational calculations, the convergence rate is good 
enough for precise extrapolation of the radiative decay 


rate. This fact is also illustrated in Fig. 10 
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